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In standard optical tomographic methods, the off-diagonal elements of a density matrix p are 
measured indirectly. Thus, the reconstruction of p, even if it is based on linear inversion, typically 
magnifies small errors in the experimental data. Recently, an optimal tomography [Phys. Rev. 

A 90 , 062123 (2014)] has been proposed theoretically to measure one-by-one all the elements of 
p. Thus, the relative errors in the reconstructed state can be the same as those in the experi¬ 
mental data. We implemented this method for two-qubit polarization states performing both local 
and global measurements. For comparison, we also experimentally implemented other well-known 
tomographic protocols based solely on local measurements (of, e.g., the Pauli operators and James- 
Kwiat-Munro-White projectors) and those with mutually unbiased bases requiring both local and 
global measurements. We reconstructed seventeen of two-qubit polarization states including sepa¬ 
rable, partially and maximally entangled. Our experiments show the highest stability against errors 
of our method in comparison to the other quantum tomographies. In particular, we demonstrate 
that each optimally-reconstructed state is embedded in the uncertainty circle of the smallest radius, 
both in terms of the trace distance and disturbance. We explain how to estimate experimentally the 
uncertainty radii for all the implemented tomographies and show that, for each reconstructed state, 
the relevant uncertainty circles intersect indicating the approximate location of the corresponding 
physical density matrix. 

PACS numbers: 03.65.Wj, 03.67.-a, 42.50.Ex 


Introduction. — Quantum tomographic methods are in¬ 
dispensable tools in experimental quantum physics. In¬ 
deed, characterizing quantum states and quantum pro¬ 
cesses are essential problems in studying the performance 
and evolution of quantum systems [T] and in develop¬ 
ing quantum technologies [5]. Moreover, both these 
problems are mathematically equivalent and are usually 
solved by applying quantum state tomography (QST). 
This approach is typically based on linear inversion [3] 
and maximum-likelihood estimation HHIl]. There are 
also other proposals of quantum state estimation based 
on, e.g., least-squares inversion [TS] as well as Bayesian 
mean estimation [Diiiiin], or linear regression estima¬ 
tion [18] . There exist dozens of QST protocols even in the 
special case of reconstructing the photonic polarization 
states (for a review see Ref. [TS] and also, e.g., Refs. pPl - 
[5Sj). Thus, choosing the best of them appeared to be 
not a simple task. However, a recent paper |29j describes 
an optimal QST protocol that minimizes the condition 
number k that characterizes the robustness against ex¬ 
perimental errors. Condition numbers were also used for 
investigating the error stability of optical tomographic 
protocols in Refs. imisniisii- In this Letter we experi¬ 
mentally study this optimal protocol in comparison to 
four other popular approaches using the same experi¬ 
mental setup for the reconstruction of two polarization- 


entangled photons. These protocols include those based 
solely on the local measurements of (i) the James-Kwiat- 
Munro-White (JKMW) projectors |1], (ii) the Pauli op¬ 
erators, and (iii) their eigenstates (the so-called stan¬ 
dard basis) Uni [31|, together with (iv) the protocol of 
Adamson and Steinberg |22j based on mutually-unbiased 
bases (MUB) and applying both local and global mea¬ 
surements, analogously to the optimal protocol. To com¬ 
pare these protocols we first derive a relation between the 
radius of an error circle associated with the reconstructed 
state and measured quantities. The radii correspond to 
the trace distance between the ideal density matrices and 
the reconstructed noisy ones. However, they can also be 
interpreted in terms of fidelity (or disturbance). 

All the approaches analyzed here are based on solving 
a linear-system problem Ax = b, where A is referred to as 
the coefficient matrix, b is the observation veetor contain¬ 
ing the measured data, and x = vec(p) is a real vector 
describing the unknown state p to be reconstructed. Here 
we choose 

x = vec(p) = [pii,Repi 2 ,Impi 2 ,Repi 3 ,Impi 3 , ...,/344]^. 

Conversely, a two-qubit density matrix p can be repre¬ 
sented as a real vector x = {xi, ...,Xiq) with its elements 
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given as follows 


p{x) 


Xi X 2 + ixs 
X2 — iX3 Xs 
X4 — ZX5 Xg — 

Xq — iX7 Xii — ixi2 


X 4 + iX 3 
Xg + ixio 
Xl3 

Xi4 "^Xi^ 


Xq + iX 7 
Xii + ixi 2 

Xi4 iXi3 

X16 


The already mentioned condition number k depends only 
on A, i.e., the system of equations used to estimate the 
density matrix from the experimental data b. The relia¬ 
bility of the reconstructed density matrix p correspond¬ 
ing to the vector x = A~^b for a given set of rotations A 
(representing our linear tomographic system) and for the 
measured data b depends on the value of k. To show the 
operational (or physical) importance of condition num¬ 
bers more explicitly, let us recall a well known theorem 
(Theorem 8.4 in Ref. [ 35 ]): Consider the system Ax = b 
with nonsingular A. Assume perturbations Sb in b. If per¬ 
turbations Sx are defined implicitly by A(a;-|-(5 x) = b+S b, 
then it holds [ 33 ] : 


k{A) ||5|| - INI - ^ ^ ||6|| ■ 


( 1 ) 


Thus, if the condition number k(A), for a given norm, is 
equal (or very close) to one, then small relative changes 
in the observation vector b imply equally small relative 
changes in the reconstructed state x. Here we calcu¬ 
late k{A) = cond 2 (A) = crmax(A)/tTniin(A) based on the 
spectral norm ||A|| 2 , which is compatible with Euclidean 
distance || • || used for other quantities in Eq. Q. The 
norm is defined by the largest singular value of A, i.e, 
IIAII 2 = max[svd(A)] = CTmax(A), where the function 
svd(A) returns the singular values of A. As shown in 
Ref. [33], the optimal tomography provides k{A) = 1 
for 16 local and nonlocal measurements. By contrast to 
this, the JKMW tomography [4] leads to k{A) = \/60.1 
for 16 local measurements, the standard separable basis 
[inii33] yields k{A) = 3 for 36 local measurements, and 
mutually-nnbiased bases [33[ [M] gives k{A) = y/E for 
20 local and nonlocal measurements. The tomography 
based on Pauli matrices gives k{A) = \/2 for 16 local 
measurements. This suggest that the density matrices 
reconstructed with these different protocols reside inside 
uncertainty circles of various radii that depend on k. 


Error analysis .— From the linearity of the linear in¬ 
version problem and Eq. Q it follows that 

1 ^ H^^ll II^^H .21 

NA) \\b + Sb\\ - N + fe|| - ^ ’\\b + SbW ^ ’ 

Let us quantify the quality of a tomography protocol with 
the trace distance E = T[p{x),p{x + fe)] = ^Tr^/JJp)^, 
where 6 p = p{6x), between the ideal, p{x), and per¬ 
turbed, p{x + Sx), density matrices. The trace distance is 
a proper measure of the distance between two density ma¬ 
trices and has a statistical interpretation in terms of the 


probability of distinguishing between the two matrices. 
Moreover, it provides a single number that quantifies the 
error introduced by the protocol. We can relate E to ||(5a;|| 
by using a standard inequality between the quadratic and 
arithmetic means of eigenvalues of (Sp)'^. The result is 

2E; = Try/Np)2 < N2d||fe||, (3) 


where, for a two-qubit density matrix, Tr[((5 p)^] = 
2 X]i=i ^xj — {5xi + 6 x 3 + 5 x \3 + feig) and the matrix di¬ 
mension is d = 4 . By combining inequalities in Eqs. (§ 
and ([^ we arrive at 


E < 



\\b + 6 b\\ 


( 4 ) 


where the random deviations Sb can be related to the 
vector of standard deviations a(b) associated with the 
mean values b. The distribution of random photon counts 
b + Sb is usually described by the Poisson statistics. Af¬ 
ter performing the measurements b+S b one assumes that 
b + S b ^ b — cr^(6), i.e., the measurement outcomes are 
the most likely (the mean) number of counts. The rela¬ 
tive error of such approximation is small if the number 
of counts is large. In order to compare the robustness of 
the tomographies, the deviations Sb need to be bounded 
from above. For the Poisson distribution ,the probability 
of a magnitude of a random deviation, \Sbi\ > 2 \/ 2 bi, 
is given by its cumulative distribution function (CDF) 
as Pr(|d6i| > 2\/Wi) = CDF(a:+) — CDF(a;_), where 
x± = [bi ± 2^/2Ei\ and for the Poisson distribution 
CDF(x < 0) = 0. The probability Pr(|(56i| > 2\/2bi) > 
0.981 is very high for all For bi > 20 its value is al¬ 
ready Pr(|d6i| > 2i/25i) > 0.993. The same approach 
applied to the Gaussian distribution provides the widely 
used 3(7 rule, which tells us that almost certainly (with 
probability 0.997) \Sbi\ < 'Sa{bi). The statistically justi¬ 
fied inequality |56i| < 2 'j 2 a{bi) leads to 


E <R = 2 VdK{A) 


||cr(b)||N + fe|| 
\\b + Sb\\ 


( 5 ) 


We have defined the uncertainty radius of the state es¬ 
timation i?, which is the maximal distance between the 
state and its estimate, in terms of only the directly mea¬ 
sured quantities. This is an important result as it allows 
to directly estimate the quality of a given state recon¬ 
struction in a very convenient way. As we will demon¬ 
strate, it also allows to visually compare the outcomes of 
various tomographies. This result allows to easily char¬ 
acterize the quality of reconstruction without knowing p 
a priori. Moreover, R can be used as a sanity check for 
the results of maximum likelihood methods, because the 
proper density matrices should be contained within the 
uncertainty circle of a radius R. 

However, using the uncertainty radius R, one can over¬ 
estimate the value of the disturbance E. Let us introduce 
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FIG. 1. Experimental setup performing both local and 
nonlocal polarization projections for the four studied tomo¬ 
graphies. Linear-optical components are quarter-wave plate 
(QWP), half-wave plate (QWP), motorized translation (MT) 
to stabilize two-photon overlap, horizontally retractable bal¬ 
anced beam splitter (BS), polarizing cube (POL). 


k = ||<56||/||cr(5)||, where 0 < fc < 2-\/2 for the Poissonian 
statistics. In the most general case we can write 


kR kR 

—F- < E< 

A^fdn^iA) ~ ~ 2V2 


( 6 ) 


where the lower bound is derived with help of Eqs. @ 
and ([^, and the relation between the trace distance 
and the Hilbert-Schmidt distance D^isip, P + ^ p) = 
•\/Tr[((5 p)2] < 2E. For an arbitrary distribution of the 
results bi, the Chebyshev’s inequality implies that the 
probability of finding the reconstructed state inside of an 
error circle of the radius kR/2^j2 is bounded from below 
by 1 — 1/A:^. This means that the minimum of 50% of 
values must lie within the y/2 standard deviations of the 
mean regardless of the distribution, i.e., the value of R/2 
bounds from above the median of E for any distribution. 
For the Poisson distribution we can find a tighter up¬ 
per bound on the probability of bi + 6 bi > bi + ka{bi) 
than the one provided by the Chebyshev inequality, i.e., 
Pr(X > x) < e~^{ep/xY (Theorem 5.4 in Ref. [35]), 
where X = bi + 5bi is the random variable, p = bi and 

X = p + k^P" 

For characterizing the quality of tomographic protocols 
we can also introduce the disturbance Db{p,P + Sp) = 
1 - [TrA/v^(/o + dp)y/p]2 = 1 - F{p,p + Sp), where 
F is the fidelity related to Bures metrics, which fulfils 
DBip,p + 6p) < T{p,p + 6p). Thus, Db < E < R. 
The disturbance was used in Ref. [22] for comparing the 
results of two-qubit tomographies. However, in our anal¬ 
ysis we apply the trace distance as it provides a more 
convenient theoretical framework. 


Experimental results .— The results obtained in the 
previous section suggest that the error E of state estima¬ 
tion depends both on k{A) and the measured quantities 
bi- Therefore, in order to compare the above-mentioned 


tomographic protocols we have prepared 17 two-qubit 
states and performed four tomographic protocols on each 
of them. These states are two-photon states described in 
the polarization basis {E[H,HV,VE[,VV}. This means 
that, e.g., pii — xi is the probability of detecting the 
two photons in the polarization state \HH) (both pho¬ 
tons are polarized horizontally). We have generated sep¬ 
arable and polarization-entangled photon pairs using the 
process of spontaneous parametric down-conversion oc¬ 
curring in a pair of BBO crystals (the so-called Kwiatet 
al. source |Mj)- With our source, we have observed about 
2.10^ two-photon detections per second having 200 mW 
of pumping power at 355 nm. Generated photons were 
subsequently brought to the input of our tomography 
setup (see Fig. [^, where the required states were pre¬ 
pared and subsequently both the local and nonlocal po¬ 
larizations projections were performed. The preparation 
of a given state was achieved using pairs of half (HWP) 
and quarter (QWP) wave plates in the input mode of 
each photon. 

In order to perform local projections on individual pho¬ 
tons, we have shifted the beam splitter BS horizontally so 
that the reflections are no longer coupled to the output 
ports. Then for each local projection, we have adjusted 
the HWP and QWP in each photon’s path and then sub¬ 
jected the photons to a polarizing cube. Respective two- 
photon detections were registered for 5 seconds. 

The nonlocal projections are achieved by combining 
the local state transformations using the HWPs and 
QWPs with the singlet-state projection on a balanced 
beam splitter BS. For this procedure to work, an addi¬ 
tional HWP (set to 45°) has to be placed in one output 
mode of the beam splitter before the photons are sub¬ 
jected to the polarizers. Again, the two-photon detec¬ 
tions were counted for 5 seconds. 

While evidently the beam splitter is superfluous for lo¬ 
cal projections and the polarizing cubes are unnecessary 
for the nonlocal projections, we maintain all the com¬ 
ponents in the setup for all times deliberately since we 
need to compare the observed detection rates across local 
and nonlocal measurements. This would be problematic 
without keeping all the components in the setup since 
the components introduce different technological losses 
(e.g. back-reflections, scattering). Further to that, our 
setup allows to switch between the local and nonlocal 
projections without much effort. 

For each tomography we have gathered the coincidence 
counts b + Sb for the specific projectors (see the Supple¬ 
ment m)- After performing the measurements we es¬ 
timated the standard deviations as b + Sb fsi b = a^b). 
This is justified for large values of bi as the relative error 
of estimating a{bi) from bi + Sbi is < {2\/2j\/bi)^^'^. In 
our experiment, we observe on average that bi « 10^ and 
the smallest values of (t( 6Q do not contribute much to 
||(T(6)|j. Thus, in order to correct for possible underes¬ 
timation of ||cr(6)||, we rescale this value by a factor of 
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FIG. 2. Experimentally recovered range of possible errors E 
for five tomographies (including the optimal, MUB, standard- 
separable-basis, Pauli operators, and JKMW protocols) for 
17 different two-qubit states. The shaded areas correspond 
to the most probable range of the error E, given by Eq. Q 
for k — y/2. The maximum error R is twice the upper limit 
r = R/2 oi the plotted error range [see Eq. B- 


1.3. In most of the tomographies, the observation vec¬ 
tor 6 -b (55 is measured directly. However, for the optimal 
tomography |29j there are 12 measurements which corre¬ 
spond to the difference of two coincidence counts, say Ci 
and c.'. In these cases, the corresponding entries of the 
observation vector 5 are bi = [{ci — c')/2j are described 
by the Skellam distribution (a generalized Poisson distri¬ 
bution), where a'^ibi) = \_{ci + c')/2j. For the Skellam 
distribution (similarly to the Poisson distribution) the 
behaviour of its cumulative distribution function implies 
that the largest disturbance can be limited (with proba¬ 
bility > 0.993 for a, c' > 20) by 55i < [2-\/2(T(5i)J. The 
results of our analysis are shown in Fig. The figure 
demonstrates that the error range E increases with the 
condition number k while using the same setup. This sug¬ 
gests that the optimal tomography can truly be the best. 
A more convincing evidence of this observation is pro¬ 
vided by analyzing the relative (trace) distances between 
the matrices reconstructed by various tomographies and 
the sizes of their error circles r = R/2. We know that 
the unperturbed density matrix must be found in the in¬ 
tersection of the error circles. This intersection is very 
close to the result of the optimal tomography because it 
has the smallest error radius. Four selected representa¬ 
tive examples of such geometric construction are shown 
in Fig. 1^ Our results for all the reconstructed states can 
be found in the Supplement m- 


FIG. 3. Relative trace distances between points correspond¬ 
ing to optimal tomography (O), standard 36 state tomog¬ 
raphy (S), MUB-based tomography (M) representing recon¬ 
structed density matrices and their corresponding disks of the 
maximum errors R for four selected reconstructed states. All 
the graphically represented distances scaled in the units of 
trace distance. The four reconstructed states can be approxi¬ 
mated with pn = \'tpn.)('ipn\, where ji/’r) = i\DR) —i\AL))/y/2, 

|V>7) = \HV), IV-9) = {\HV) - \VH))/V2, |V>i4) = lei.eii,), 
where jeia) = (-0.6556-1-0.62481)117)-1-0.4241 |F) and jei;,) = 
(-0.1415-0.71651)177)-b0.6831|F), where |77), \V), \D), |A), 

17?), and |L) stand for the horizontal, vertical, diagonal, an¬ 
tidiagonal, left-circular, and right circular states, respectively. 
An ideally reconstructed state lies in the intersection of all the 
error disks of radius R. Note that the discs of radius R/2, as¬ 
sociated with the most probable range of errors E (see Fig.j^, 
do not necessarily intersect. 

Conclusions. — We have for the first time implemented 
the optimal two-qubit tomography and compared it with 
the other four important tomographic protocols. This 
method corresponds to measuring one by one all the el¬ 
ements pnm of a density matrix p. This is in contrast 
to the other protocols, where the off-diagonal elements 
of p are measured indirectly, i.e., the measured photon 
numbers correspond to linear combinations of some ele¬ 
ments Pnm- We have developed a method for estimating 
the error radii (in units of the trace distance) of circles 
containing the reconstructed density matrices. We have 
demonstrated that all the tomographies can be imple¬ 
mented and compared using the same framework. Our 
results confirm that the optimal tomography provides the 
most reliable results among all other analyzed protocols. 
This makes the optimal tomography a method of choice 
if the quality of the reconstructed density matrix is a 
priority. 
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Priority Choice Experimental Two-qubit Tomography: 

Measuring One by One All Elements of Density Matrices: 
Supplementary Materials 

Here we show explicitly all the density matrices discussed in the Letter, which are reconstructed with the optimal 
tomographic protocol and those based on: (i) mutually unbiased bases, (ii) the James-Kwiat-Munro-White projec¬ 
tors, (iii) the tensor products of the Pauli operators, and (iv) the standard separable basis corresponding to all the 
eigenvectors of the Pauli operators. We also present the coefhcient matrices, observation vectors corresponding to 
coincidence counts, the estimated variances for the observations, and the error radii for each reconstructed matrix. 
Finally, we compare the reconstructed matrices graphically, where we show the relative trace distances between the 
reconstructed states and they error radii. 

Reconstructed density matrices 

The 17 density matrices are reconstructed by solving linear inversion problem for four tomographies. We have 
prepared 17 different states of high purity, which approximately correspond to: 

lyp = {\HH) -\VV))/V2, \4,2) = {\HH) + \VV))/V2, = i\HH) - i\VV))/V2, 

\^l,i) = {\DR)-^\AL))/V2, \^,) = {\HV)+i\VH))/V2, \i,^) = {\HV) A\VH))/V2, 

\^Pj) = \HV), \^^) = {\HH)+i\VV))/^/2, lyg) = (lidP) - |yi7))/V2, 

l^io) = {\HV) - i\VH))/V2, |V>ii) = {\DL) + i\AR))lV2, 1 ^ 12 ) = {\DL) - i\AR))/V2, 

Ibis) = leiaeifc), Ibi4) = \e 2 ae 2 b), Ibis) = 0.79|i7P) - 0.61|PiJ), 

|,/,ie) = 0.50|i7P) - 0.87|Pi7), Ibir) = 0.35|i7y) - 0.94|Pi7); 

where |eia) = (-0.6556 -f 0.6248b|i7) -h 0.4241|P), \eib) = (-0.1415 - 0.7165b|i7) -h 0.6831|P), psa) = (-0.9608 -h 
0.2091i)|i7)-|-0.1822|P), and |e 2 h) = (0.2613-1-0.7338i)|i7)-|-0.6271|P) are single photon elliptic polarization states. We 
mark the data relevant to a particular tomography as follows: index O for the optimal tomography; S for the standard 
36-state tomography; J for the James-Kwiat-Munro-White (JKMW) protocol; M for the MUB-based tomography; 
P for the Pauli matrices based tomography. 

Standard 36-state tomography 


0.4922 0.0020+ 0.0156i 

0.0020 - 0.0156Z 0.0047 

-0.0042 - 0.0354t -0.0054 - 0.0228z 
-0.4607-k 0.0750Z 0.0255 -f 0.0002i 


-0.0042 -t 0.0354i -0.4607 - 0.0750t 
-0.0054 -t 0.0228Z 0.0255 - 0.0002t 
0.0136 0.0184-k 0.0656Z 

0.0184 - 0.0656Z 0.4895 


0.4870 0.0029 - 0.0358z 0.0237 - 0.0103z 

0.0029 -f 0.0358Z 0.0085 -0.0000 - 0.0219z 

0.0237-f 0.0103Z -0.0000-f 0.0219t 0.0038 

0.4723 - 0.0515Z 0.0244 - 0.0413z -0.0052 -k 0.0378z 


0.4723-k 0.0515Z 
0.0244-k 0.0413Z 
-0.0052 - 0.0378t 
0.5007 


PS,3 = 


0.5363 0.0744-O.Oeilz 0.0830 - 0.0243z -0.0027-k 0.4636z 

0.0744-k 0.0611Z 0.0206 0.0522 - 0.0036z -0.0513-k 0.0602t 

0.0830-k 0.0243Z 0.0522 -t 0.0036z 0.0002 0.0005-k 0.0402t 

-0.0027- 0.4636Z -0.0513 - 0.0602t 0.0005 - 0.0402t 0.4429 



PS ,4 = 


0.3005 0.2633 -f 0.0250t 0.0456 - 0.2223z 

0.2633 - 0.0250Z 0.2482 0.0417 - 0.184U 

0.0456-f 0.2223Z 0.0417-f 0.1841i 0.1270 

-0.0648 - 0.2710Z -0.0334 - 0.2402t -0.2145 -f 0.0463z 


-0.0648 -f 0.2710t 
-0.0334 -k 0.2402t 
-0.2145 - 0.0463Z 
0.3244 



PS,b = 


PSfi = 


PS,7 = 


PS,S = 


PS,9 = 


PS,10 


PS,11 = 


PS,12 


PS,13 


PS,14 = 


0.0135 

-0.0859+ 0.0013Z 

0.0334 + 0.0685Z 

0.0159+ 0.0018Z ■ 

-0.0859- 0.00131 

0.5111 

0.0397 - 0.4647Z 

-0.0002 - 0.0441i 

0.0334 - 0.0685Z 

0.0397 + 0.4647Z 

0.4697 

0.0078 - 0.0236Z 

0.0159- 0.0018Z 

-0.0002 + 0.044H 

0.0078 + 0.0236Z 

0.0056 

0.0146 

0.0450 + 0.0866Z 

0.0849 + 0.0819Z 

-0.0046 + 0.0341i' 

0.0450 - 0.0866Z 

0.4606 

0.4608 - 0.0160Z 

-0.0222 - 0.0282i 

0.0849 - 0.0819Z 

0.4608 + 0.0160Z 

0.5177 

-0.0526 - 0.0125Z 

-0.0046 - 0.0341i 

-0.0222 + 0.0282i 

-0.0526+ 0.0125Z 

0.0072 

0.0005 

0.0578 - 0.0482i 

0.0064 + 0.0028Z 

-0.0088 - 0.0005z' 

0.0578 + 0.0482i 

0.9915 

0.0290 + 0.0678Z 

-0.0436-0.1064i 

0.0064 - 0.0028Z 

0.0290 - 0.0678Z 

0.0049 

-0.0043 - 0.0022i 

-0.0088 + 0.0005Z 

-0.0436+ 0.1064i 

-0.0043 + 0.0022i 

0.0032 

0.5609 

-0.0543+ 0.0315Z 

0.0357 - 0.0364i 

0.0027 - 0.4704i ‘ 

-0.0543 - 0.0315Z 

0.0067 

-0.0470 + 0.0079Z 

-0.0016 + 0.0511Z 

0.0357+ 0.0364i 

-0.0470 - 0.0079Z 

-0.0091 

0.0065 - 0.0126Z 

_ 0.0027+ 0.4704i 

-0.0016- 0.0511Z 

0.0065 + 0.0126Z 

0.4416 

-0.0208 

-0.0024 - 0.0670Z 

-0.0036 + 0.0319i 

-0.0224 - 0.0403z' 

-0.0024 + 0.0670Z 

0.5767 

-0.4584 - 0.0718Z 

-0.0076 - 0.0478Z 

-0.0036- 0.0319Z 

-0.4584 + 0.0718Z 

0.4334 

0.0134+ 0.0045Z 

-0.0224 + 0.0403Z 

-0.0076 + 0.0478Z 

0.0134 - 0.0045Z 

0.0107 

-0.0118 

0.0243 - 0.0103Z 

0.0134 + 0.0063Z 

-0.0090 - 0.0064i' 

0.0243 + 0.0103Z 

0.5080 

0.0499 + 0.4684i 

-0.0174- 0.0050Z 

0.0134 - 0.0063Z 

0.0499 - 0.4684i 

0.4801 

0.0302 + 0.0537Z 

-0.0090 + 0.0064i 

-0.0174+ 0.0050Z 

0.0302 - 0.0537Z 

0.0237 

0.2826 

0.2502+ 0.0153Z 

-0.0157+ 0.2358Z 

0.0053 - 0.2593Z ' 

0.2502 - 0.0153Z 

0.2221 

-0.0295 + 0.2241i 

0.0222 - 0.2349Z 

-0.0157- 0.2358Z 

-0.0295 - 0.2241i 

0.2718 

-0.2432 + 0.0003Z 

0.0053+ 0.2593Z 

0.0222 + 0.2349Z 

-0.2432 - O.OOOSz 

0.2235 

0.2100 

-0.2611 - 0.0385Z 

-0.0419- 0.2516? 

0.0345 - 0.2156z' 

-0.2611 + 0.0385Z 

0.2815 

0.0452 + 0.2460Z 

0.0115 + 0.2405Z 

-0.0419 + 0.2516Z 

0.0452 - 0.2460Z 

0.2312 

0.2259 + 0.0926Z 

0.0345 + 0.2156Z 

0.0115 - 0.2405Z 

0.2259 - 0.0926Z 

0.2772 

0.2849 

-0.2477 - 0.0400Z 

0.0073 + 0.2296Z 

0.0167+ 0.2654i ' 

-0.2477 + 0.0400Z 

0.2042 

0.0104- 0.2183Z 

-0.0427- 0.2274i 

0.0073 - 0.2296Z 

0.0104+ 0.2183Z 

0.2687 

0.2297 - 0.0240i 

0.0167 - 0.2654i 

-0.0427+ 0.2274i 

0.2297 + 0.0240i 

0.2422 

0.3786 

-0.0779 - 0.3764i 

-0.1289+ 0.1049Z 

0.1243+ 0.0775Z ' 

-0.0779 + 0.3764i 

0.3745 

-0.0776- 0.1422i 

-0.0950 +O.lOOOz 

-0.1289 - 0.1049Z 

-0.0776+ 0.1422i 

0.1345 

-0.0377- 0.1163Z 

0.1243- 0.0775Z 

-0.0950- O.lOOOz 

-0.0377+ 0.1163Z 

0.1124 



PS, 15 = 


PS, 16 — 


PS,17 = 


PJ,1 = 


PJ,2 


PJ,3 


pja = 


PJ,5 = 


PJfi — 


0.5241 

0.1701 + 0.40987 

-0.0732 + 0.04977 

-0.0613 - 0.07407 

0.1701 - 0.40987 

0.3733 

0.0048 + 0.08007 

-0.0776 + 0.01707 

-0.0732 - 0.04977 

0.0048 - 0.08007 

0.0586 

0.0165 + 0.05197 

-0.0613 + 0.07407 

-0.0776- 0.01767 

0.0165 - 0.0519i 

0.0440 

0.0079 

-0.0088- 0.11917 

-0.0026 + 0.04697 

-0.0179 - 0.03697 

-0.0088 + 0.11917 

0.6443 

-0.4307- 0.05977 

-0.0144- 0.08267 

-0.0026 - 0.04697 

-0.4307 + 0.05977 

0.3254 

0.0248 + 0.02987 

-0.0179 + 0.03697 

-0.0144 + 0.08267 

0.0248 - 0.02987 

0.0224 

0.0071 

-0.0098- 0.1224i 

-0.0079 + 0.04507 

-0.0180- 0.03317' 

-0.0098+ 0.1224i 

0.7482 

-0.3769 - 0.0542i 

-0.0215 - 0.10087 

-0.0079 - 0.04507 

-0.3769 + 0.0542i 

0.2232 

0.0274 + 0.01977 

-0.0180 + 0.03317 

-0.0215 + 0.10087 

0.0274 - 0.01977 

0.0215 


JKMW Tomography 


0.4879 

-0.0241 - 0.0194i 
-0.0198 - 0.0473i 
-0.4503 + 0.0438i 


-0.0241 + 0.0194i 
0.0054 

-0.0313-0.11937 
0.0428 + 0.00607 


-0.0198 + 0.04737 
-0.0313 + 0.11937 
0.0225 

-0.0023 - 0.0852i 


-0.4503 - 0.04387 
0.0428 - 0.00667 
-0.0023 + 0.0852i 
0.4842 


0.4748 0.0156 - 0.0702i 

0.0156 + 0.0702i 0.0107 

-0.0009 - 0.00897 0.0457 + 0.1052i 
0.4543 - 0.0190i 0.0248 - 0.0104i 


-0.0009 + 0.0089i 
0.0457-0.1052i 
0.0079 

0.0098 + 0.03987 


0.4543 + 0.01907 
0.0248 + 0.0104i 
0.0098 - 0.03987 
0.5065 


0.5221 

0.0753 + 0.07077 
0.0706 + 0.02387 
0.0314- 0.44517 


0.0753 - 0.07077 
0.0213 

0.0016 + 0.03387 
-0.0532 - 0.04407 


0.0706 - 0.02387 0.0314 + 0.44517 
0.0016 - 0.03387 - 0.0532 + 0.04467 
0.0112 0.0008 + 0.06187 

0.0008 - 0.06187 0.4454 


0.2807 

0.2535 - 0.0302i 
0.0392 + 0.18537 
-0.0417- 0.25157 


0.2535 + 0.0302i 
0.2543 

-0.0142 + 0.14737 
-0.0489 - 0.19997 


0.0392 - 0.18537 
-0.0142-0.14737 
0.1256 

-0.2106 + 0.00757 


-0.0417 + 0.25157 
-0.0489 + 0.19997 
-0.2106 - 0.00757 
0.3394 


0.0260 

-0.0955 - 0.0352i 
0.0370 - 0.04507 
0.0842 - 0.01437 


-0.0955 + 0.0352i 
0.5035 

0.0139+ 0.4182i 
-0.0177+ 0.0134i 


0.0370 + 0.04507 
0.0139 - 0.4182i 
0.4676 

-0.0199 + 0.02157 


0.0842 + 0.01437 
-0.0177-0.0134i 
-0.0199 - 0.02157 
0.0029 


0.0367 

0.0389 - 0.10597 
0.0675 - 0.03957 
0.0312 - 0.01407 


0.0389 + 0.10597 0.0675 + 0.03957 

0.4500 0.4370 - 0.07107 

0.4370 + 0.07167 0.5064 

-0.0533 + 0.00997 -0.0482 + O.OMli 


0.0312 + 0.01467 
-0.0533 - 0.00997 
-0.0482 - O.OMli 
0.0069 



PJ,7 — 


PJ,S = 


PJ,9 = 


PJ,10 


Pj,11 


PJ,12 = 


PJ,13, = 


pjm = 


PJ,15 — 


Pj,16 = 


0.0062 0.0456 - 0.0314i 

0.0456 + 0.0314i 0.9818 

0.0057 - 0.0030i 0.0226 - 0.0910i 
0.0072 + 0.0079i -0.0969 + 0.0714i 

0.5550 -0.0784 + 0.0200i 

-0.0784 - 0.0200* 0.0135 

-0.0058 + 0.0372i 0.0219 - 0.0377* 
-0.0223 + 0.4409* -0.0032 - 0.0374i 

0.0089 -0.0324 - 0.0301* 

-0.0324 + 0.0301* 0.5684 

-0.0031 + 0.0189* -0.4348 - 0.0409* 
0.0009 + 0.0043* -0.0166 + 0.0446* 

0.0058 0.0400 - 0.0019i 

0.0400 + 0.0019* 0.5133 

0.0128 + 0.0240* 0.0953 - 0.3781* 

-0.0401 + 0.0454i -0.0293 - 0.0144 


0.0057 + 0.0030* 0.0072 - 0.0079* 

0.0226 + 0.0910* -0.0969- 0.0714i 
0.0032 -0.0035 + 0.0032i 

-0.0035 - 0.0032i 0.0087 

-0.0058 - 0.0372i -0.0223 - 0.4409*' 
0.0219 + 0.0377* -0.0032 + 0.0374i 
0.0043 -0.0027 - 0.0162i 

-0.0027+ 0.0162i 0.4271 

-0.0031 - 0.0189* 0.0009 - 0.0043* ' 
-0.4348 + 0.0409i -0.0166 - 0.0446* 
0.4209 0.0002- 0.0151* 

0.0002 + 0.0151* 0.0018 

0.0128 - 0.0240* -0.0401 - 0.0454i' 
0.0953 + 0.3781* -0.0293 + 0.0144i 
0.4753 0.0619 + 0.0094i 

0.0619 - 0.0094i 0.0056 


0.2918 0.2453- 0.0200* -0.0285 + 0.1767* -0.0527- 0.2169* 

0.2453 + 0.0200* 0.1992 0.0202 + 0.2358* 0.0220 - 0.2014i 

-0.0285 - 0.1767* 0.0202 - 0.2358* 0.2948 -0.2464 + 0.0008* 

-0.0527 + 0.2169* 0.0220 + 0.2014i -0.2464-0.0008* 0.2142 


0.1923 

-0.2367- 0.0390* 

-0.0338-0.1873* 

-0.0453 - 0.2798* 

-0.2367 + 0.0390* 

0.3002 

0.0996+ 0.1912i 

0.0039 + 0.2219* 

-0.0338 + 0.1873* 

0.0996- 0.1912i 

0.2127 

0.2404 + 0.0740* 

-0.0453 + 0.2798* 

0.0039 - 0.2219i 

0.2404 - 0.0740* 

0.2947 

0.3064 

-0.2350- 0.0157* 

0.0162+ 0.1740i 

0.0769 + 0.2162i ' 

-0.2350 + 0.0157* 

0.1884 

-0.0753 - 0.2582i 

-0.0327-0.1833* 

0.0162 - 0.1740* 

-0.0753 + 0.2582i 

0.2840 

0.2538 - 0.0131* 

0.0769 - 0.2162i 

-0.0327 + 0.1833* 

0.2538 + 0.0131* 

0.2211 

0.3803 

-0.0676 - 0.4027* 

-0.1333+ 0.1142i 

0.1197 + 0.1053* ■ 

-0.0676 + 0.4027* 

0.3748 

-0.0745 - 0.1464i 

-0.0932 + 0.1079* 

-0.1333- 0.1142i 

-0.0745 + 0.1464i 

0.1345 

-0.0335 - 0.121H 

0.1197- 0.1053* 

-0.0932 - 0.1079i 

-0.0335 + 0.1211i 

0.1104 

0.5276 

0.1690 + 0.4181* 

-0.0627+ 0.0920i 

-0.0903 - 0.0653* 

0.1690- 0.4181* 

0.3819 

0.0342 + 0.0807* 

-0.0764 + 0.0379* 

-0.0627- 0.0920* 

0.0342 - 0.0807* 

0.0478 

0.0032 + 0.0440* 

-0.0903 + 0.0653* 

-0.0764 - 0.0379i 

0.0032 - 0.0440i 

0.0427 

0.0207 

-0.0045 - 0.1146* 

-0.0112 + 0.0073* 

-0.0325 - 0.0020* 

-0.0045 + 0.1146* 

0.6312 

-0.4297 + 0.0089* 

0.0034 - 0.0464i 

-0.0112 - 0.0073* 

-0.4297 - 0.0089i 

0.3450 

0.0040 + 0.0183* 

-0.0325 + 0.0020* 

0.0034 + 0.0464i 

0.0040 - 0.0183* 

0.0031 

0.0202 

-0.0172 - 0.1220* 

-0.0095 + 0.0268* 

-0.0629 - 0.0175*' 

-0.0172 + 0.1220* 

0.7415 

-0.3322 - 0.0029* 

-0.0138 - 0.0503* 

-0.0095 - 0.0268* 

-0.3322 + 0.0029i 

0.2341 

0.0108 - 0.0009* 

-0.0629 + 0.0175* 

-0.0138 + 0.0503* 

0.0108 + 0.0009i 

0.0042 


Pj,17 = 



MUB-based tomography 


Pm,i = 


Pm,2 = 


Pm,3 = 


Pm,4 — 


Pm, 5 


Pm,6 = 


Pm,7 = 


Pm,8 = 


Pm,9 — 


0.4789 0.0857 + 0.02751 

0.0857-0.02751 0.0311 

0.0044 - 0.0534i -0.0110 + 0.0029i 
-0.4699 + 0.04767 0.0332 + 0.01967 


0.0044 + 0.0534i 
-0.0110- 0.00297 
-0.0015 

0.0215 - 0.07587 


0.4888 -0.0083 - 0.0456? 

-0.0083 + 0.04567 0.0296 

0.0198 - 0.00137 0.0027 + 0.0241i 

0.4919 - 0.0541i 0.0206- 0.02917 


0.0198 + 0.00137 
0.0027- 0.02417 
-0.0322 

-0.0820 + 0.0477? 


0.5214 

0.0980 + 0.0452i 
0.0865 + 0.0535? 
0.0541 - 0.4566? 


0.0980 - 0.0452i 
0.0573 

-0.0021 - 0.0016? 
-0.0472 - 0.0870? 


0.0865 - 0.0535? 
-0.0021 + 0.0016? 
0.0053 

-0.0563 - 0.0557? 


0.2423 

0.2560 - 0.0449? 
0.0503 + 0.2014i 
0.0328 - 0.2533? 


0.2560 + 0.0449? 
0.3304 

0.1157 + 0.1609? 
-0.0336- 0.2319? 


0.0503 - 0.2014i 
0.1157- 0.1609i 
0.1676 

-0.1754 + 0.0308? 


-0.0004 

-0.0941 - 0.0188? 
0.0384 - 0.0899i 
0.0017 + 0.0313? 


-0.0941 + 0.0188? 
0.4975 

0.1830+ 0.4264i 
0.0052 + 0.0638? 


0.0384 + 0.0899? 
0.1830- 0.4264i 
0.5160 

0.0302 + 0.0058? 


0.0119 

0.0116 - 0.0964i 
0.0767 - 0.0744i 
0.0106 - 0.0150? 


0.0116+ 0.0964i 
0.4477 

0.4659- 0.0031? 
-0.0301 + 0.0330? 


0.0767 + 0.0744i 
0.4659 + 0.0031? 
0.5525 

-0.0780 + 0.0024i 


0.0377 

0.0294 + 0.0443? 
0.0050 + 0.0081? 
-0.0075 - 0.0315? 


0.0294 - 0.0443? 
0.9628 

0.0006 - 0.0373? 
-0.0462+ 0.1064i 


0.0050- 0.0081? 
0.0006 + 0.0373? 
0.0047 

0.0296 - 0.0029? 


0.5549 

-0.0058 - 0.0214i 
0.0324 - 0.0106? 
-0.0598 + 0.4494i 


-0.0058 + 0.0214i 
0.0236 

0.0055 + 0.0171? 
-0.0052 - 0.0074i 


0.0324 + 0.0106? 
0.0055 - 0.0171? 
-0.0329 

-0.0199 + 0.0231? 


-0.0157 

-0.0020 + 0.0764i 
-0.0028 - 0.0180i 
-0.0136+ 0.0240i 


-0.0020 - 0.0764i 
0.5405 

-0.4696 + 0.0885? 
-0.0067 + 0.0289? 


-0.0028 + 0.0180? 
-0.4696 - 0.0885? 
0.4525 

0.0798 + 0.0047? 


-0.4699 - 0.0476? 
0.0332 - 0.0196? 
0.0215 + 0.0758? 
0.4915 


0.4919+ 0.0541i 
0.0206 + 0.0291? 
-0.0820 - 0.0477? 
0.5138 


0.0541 + 0.4566? 
-0.0472 + 0.0870? 
-0.0563 + 0.0557? 
0.4160 


0.0328 + 0.2533? 
-0.0336 + 0.2319? 
-0.1754- 0.0308? 
0.2597 


0.0017- 0.0313? 
0.0052 - 0.0638? 
0.0302 - 0.0058? 
-0.0131 


0.0106 + 0.0150? 
-0.0301 - 0.0330? 
-0.0780 - 0.0024i 
- 0.0121 


-0.0075 + 0.0315? 
-0.0462 - 0.1064i 
0.0296 + 0.0029? 
-0.0052 


-0.0598 - 0.4494i 
-0.0052 + 0.0074i 
-0.0199- 0.0231? 
0.4544 


-0.0136 - 0.0240? 
-0.0067- 0.0289? 
0.0798 - 0.0047? 
0.0228 




-0.0061 

-0.0415 - 0.0233Z 

0.0029 - 0.0435Z 

-0.0009 - 0.0034Z 

Pm,10 = 

-0.0415 + 0.0233i 

0.4724 

-0.0159 + 0.4682i 

-0.0280 + 0.0304i 

0.0029 + 0.0435Z 

-0.0159 - 0.4682i 

0.4984 

0.0302 + 0.0410Z 


-0.0009 + 0.0034i 

-0.0280 - 0.0304i 

0.0302 - 0.0410Z 

0.0353 


0.3090 

0.2536 - 0.0020Z 

-0.0236 + 0.2387Z 

0.0049 - 0.2545Z ' 

Pm,11 = 

0.2536 + 0.0020i 
-0.0236 - 0.2387i 

0.2169 

-0.0266 - 0.2206Z 

-0.0266 + 0.2206Z 

0.2625 

0.0128 - 0.2344i 
-0.2126 - 0.0164i 


0.0049 + 0.2545Z 

0.0128 + 0.2344i 

-0.2126 + 0.0164i 

0.2116 


0.2210 

-0.2802 - 0.0510Z 

-0.0475 - 0.2454i 

-0.0166 - 0.2218z' 

Pm,12 = 

-0.2802 + 0.0510Z 

0.2512 

0.0112 + 0.2530Z 

0.0076 + 0.2359Z 

-0.0475 + 0.2454i 

0.0112 - 0.2530Z 

0.1993 

0.1977 + 0.0841Z 


-0.0166 + 0.2218i 

0.0076 - 0.2359Z 

0.1977-0.0841i 

0.3285 


0.3059 

-0.2612 - 0.0147Z 

0.0146 + 0.2246Z 

0.0422 + 0.2544i ' 

Pm,13 = 

-0.2612 + 0.0147i 
0.0146 - 0.2246Z 

0.1762 

-0.0175 + 0.2085Z 

-0.0175 - 0.2085Z 
0.2792 

-0.0342 - 0.2370Z 
0.1911 + 0.0008Z 


0.0422 - 0.2544i 

-0.0342 + 0.2370Z 

0.1911 -0.0008Z 

0.2387 


0.4095 

-0.0753 - 0.3984i 

-0.1351 + 0.1251Z 

0.1159 + 0.0932i ■ 

PM,14 = 

-0.0753 + 0.3984i 

0.3649 

-0.0998-0.1619i 

-0.0991 + 0.1004i 

-0.1351 - 0.1251Z 

-0.0998+ 0.1619Z 

0.1005 

-0.0384-0.1226Z 


0.1159 - 0.0932i 

-0.0991 - 0.1004i 

-0.0384 + 0.1226Z 

0.1251 


0.5590 

0.1688 + 0.4115Z 

-0.0750 + 0.0461Z 

-0.0342 - 0.0741Z 

PM,15 = 

0.1688-0.4115Z 

0.3472 

0.0472 + 0.0803Z 

-0.0795 + 0.0198Z 

-0.0750 - 0.0461i 

0.0472 - 0.0803Z 

0.0258 

0.0397 + 0.0409Z 


-0.0342 + 0.0741i 

-0.0795 - 0.0198Z 

0.0397 - 0.0409Z 

0.0680 


0.0095 

-0.0057-0.1331Z 

-0.0030 + 0.0330i 

-0.0239 - 0.0074Z 

PM,16 = 

-0.0057 + 0.1331i 
-0.0030 - 0.0330Z 

0.6278 

-0.3716 + 0.0931Z 

-0.3716 - 0.0931Z 
0.3255 

-0.0152 - 0.0669Z 
0.0919 + 0.0219Z 


-0.0239 + 0.0074i 

-0.0152 + 0.0669Z 

0.0919 - 0.0219Z 

0.0372 


0.0219 

-0.0122 - 0.1428Z 

-0.0034 + 0.0132i 

-0.0230 - 0.0034i 

PM,17 = 

-0.0122 + 0.1428i 

0.7250 

-0.3510 - 0.0873Z 

-0.0176 - 0.0751Z 

-0.0034- 0.0132i 

-0.3510 + 0.0873Z 

0.2151 

0.0860 + 0.0050Z 


-0.0230 + 0.0034i 

-0.0176 + 0.0751Z 

0.0860 - 0.0050Z 

0.0379 


Optimal tomography 


0.4879 -0.0194 + 0.0278i 

-0.0194 - 0.0278Z 0.0054 

0.0045 - 0.0413Z -0.0112 + O.OOOSz 
-0.4760 + 0.0086Z 0.0336 - 0.0064i 


0.0045 + 0.0413Z -0.4760 - 0.0086z 
-0.0112 - 0.0003i 0.0336 + 0.0064i 
0.0225 -0.0033 + 0.0768z 

-0.0033 - 0.0768i 0.4842 


Po,i = 



PO,2 


PO,3 = 


POA = 


PO,5 = 


POfi = 


PO,7 = 


POA = 


PO,9 — 


PO,10 = 


0.4748 

0.0191 + 0.04431 
0.0193 + 0.01301 
0.4781 - 0.00191 


0.0191 - 0.0443i 
0.0107 

0.0026- O.OOOH 
0.0200 - 0.0379* 


0.0193 - O.OlSOi 0.4781 + 0.0019* 
0.0026 + 0.0001* 0.0200 + 0.0379* 
0.0079 0.0111 - 0.0464i 

0.0111+ 0.0464i 0.5065 


0.5221 0.0816 - 0.0442i 

0.0816 + 0.0442i 0.0213 

0.0847 + 0.0240* -0.0020 + 0.0026* 
0.0530 - 0.4532i -0.0463 - 0.0584i 


0.0847 - 0.0240i 
-0.0020 - 0.0026* 
0.0112 

0.0095 - 0.0545* 


0.0530 + 0.4532i 
-0.0463 + 0.0584i 
0.0095 + 0.0545* 
0.4454 


0.2807 0.2569 + 0.0406* 

0.2569 - 0.0406* 0.2543 

0.0456 + 0.2118* 0.1047 + 0.1856* 
0.0297- 0.2155* -0.0304- 0.2331* 


0.0456 - 0.2118* 
0.1047-0.1856* 
0.1256 

-0.2027 + 0.0279i 


0.0297 + 0.2155* 
-0.0304 + 0.2331* 
-0.2027 - 0.0279* 
0.3394 


0.0260 -0.0969 + 0.0186* 

-0.0969 - 0.0186* 0.5035 

0.0380- 0.0730* 0.1813 + 0.4222i 

0.0017 + 0.0078* 0.0051 + 0.0374i 


0.0380 + 0.0730i 0.0017 - 0.0078* 

0.1813 - 0.4222i 0.0051 - 0.0374i 

0.4676 -0.0051 - 0.0057* 

-0.0051 + 0.0057* 0.0029 


0.0367 0.0502 + 0.0923* 

0.0502 - 0.0923* 0.4500 

0.0734 - 0.0820* 0.4460 + 0.0748* 
0.0102 + 0.0178* -0.0288 + 0.0230i 


0.0734 + 0.0820* 0.0102-0.0178* 

0.4460 - 0.0748* -0.0288 - 0.0230* 
0.5064 -0.0430 - 0.0023* 

-0.0430 + 0.0023* 0.0069 


0.0062 

0.0557+ 0.0422i 
0.0047 - 0.0036* 
-0.0072 - 0.0029* 


0.0557- 0.0422i 
0.9818 

0.0006 - 0.0380i 
-0.0441 + 0.1031* 


0.0047 + 0.0036* 
0.0006 + 0.0380i 
0.0032 

-0.0049 - 0.0027* 


-0.0072 + 0.0029* 
-0.0441 -0.1031* 
-0.0049 + 0.0027* 
0.0087 


0.5550 

-0.0622 - 0.0204i 
0.0308 + 0.0306* 
-0.0569 + 0.4699* 


-0.0622 + 0.0204i 
0.0135 

0.0053 + 0.0015* 
-0.0050 - 0.0533* 


0.0308 - 0.0306* 
0.0053 - 0.0015* 
0.0043 

-0.0039 + 0.0220* 


-0.0569 - 0.4699* 
-0.0050 + 0.0533* 
-0.0039 - 0.0220* 
0.4271 


0.0089 

-0.0219 + 0.0731* 
-0.0027- 0.0214i 
-0.0130 - 0.0040* 


-0.0219- 0.0731* 
0.5684 

-0.4492 - 0.0011* 
-0.0065 + 0.0551* 


-0.0027+ 0.0214i 
-0.4492 + 0.0011* 
0.4209 

-0.0068 + 0.0044i 


-0.0130 + 0.0040* 
-0.0065 - 0.0551* 
-0.0068 - 0.0044i 
0.0018 


0.0058 0.0406 - 0.0232i 

0.0406 + 0.0232i 0.5133 

0.0029 - 0.0025* -0.0158 - 0.4665* 
-0.0009 - 0.0092i -0.0279 + 0.0088* 


0.0029 + 0.0025* -0.0009 + 0.0092i 
-0.0158 + 0.4665* -0.0279 - 0.0088* 
0.4753 0.0465 + 0.0408* 

0.0465 - 0.0408* 0.0056 


0.2918 0.2395 - 0.0021i -0.0242 + 0.2258* 

0.2395 + 0.0021i 0.1992 -0.0274 + 0.2259* 

-0.0242 - 0.2258* -0.0274- 0.2259* 0.2948 

0.0050 + 0.2663* 0.0132 + 0.2382i -0.2469 + 0.0169* 


0.0050 - 0.2663* 
0.0132 - 0.2382i 
-0.2469 - 0.0169* 
0.2142 


Po,ii = 




0.1923 

-0.2541 - 0.0502i 

-0.0467 - 0.24701 

PO,12 = 

-0.2541 + 0.0502i 
-0.0467 + 0.2470* 

0.3002 

0.0111 - 0.2257* 

0.0111 + 0.22571 
0.2127 


-0.0163 + 0.2401* 

0.0075 - 0.2518* 

0.2396 - 0.08271 


0.3064 

-0.2359- 0.0153* 

0.0152 + 0.22701 

PO,13 = 

-0.2359 + 0.0153* 

0.1884 

-0.0182- 0.19921 

0.0152 - 0.2270* 

-0.0182 + 0.1992i 

0.2840 


0.0439- 0.2918* 

-0.0356 + 0.2367* 

0.2487 - 0.00091 


0.3803 

-0.0814 - 0.3887* 

-0.1318 + 0.10511 

POM = 

-0.0814 + 0.3887* 
-0.1318-0.1051* 

0.3748 

-0.0973 + 0.1265* 

-0.0973- 0.12651 
0.1345 


0.1131 - 0.0665* 

-0.0967- 0.1000* 

-0.0398 + 0.11961 


0.5276 

0.1680 + 0.4105* 

-0.0748 + 0.05271 

POM = 

0.1680- 0.4105* 

0.3819 

0.0471 + 0.08741 

-0.0748 - 0.0527* 

0.0471 - 0.08741 

0.0478 


-0.0341 + 0.0538* 

-0.0793- 0.0196* 

0.0094 - 0.04081 


0.0207 

-0.0261 - 0.1296* 

-0.0029 + 0.04011 

POM = 

-0.0261 + 0.1296* 

0.6312 

-0.3619- 0.00961 

-0.0029 - 0.0401* 

-0.3619 + 0.0096* 

0.3450 


-0.0233 + 0.0012i 

-0.0148 + 0.09121 

0.0080 - 0.02131 


0.0202 

-0.0269- 0.13851 

-0.0033 + 0.04861 

POM = 

-0.0269 + 0.1385* 
-0.0033 - 0.0486* 

0.7415 

-0.3405 - 0.01021 

-0.3405 + 0.01021 
0.2341 


-0.0223 - 0.0155* 

-0.0171 + 0.09851 

0.0106 - 0.00491 


Pauli matrices based tomography 


0.4879 

-0.0194 + 0.02781 

0.0045 + 0.04131 

-0.0194 - 0.02781 

0.0054 

-0.0053 + 0.02231 

0.0045 - 0.04131 

-0.0053 - 0.02231 

0.0225 

-0.4516 + 0.07351 

0.0336 - 0.00641 

-0.0033 - 0.07681 

0.4748 

0.0191 - 0.04431 

0.0193 - 0.01301 

0.0191 + 0.04431 

0.0107 

-0.0000 - 0.02161 

0.0193 + 0.01301 

-0.0000 + 0.02161 

0.0079 

0.4662 - 0.05081 

0.0200 - 0.03791 

0.0111 + 0.04641 

0.5221 

0.0816 - 0.04421 

0.0847 - 0.02401 

0.0816 + 0.04421 

0.0213 

0.0509 - 0.00351 

0.0847 + 0.02401 

0.0509 + 0.00351 

0.0112 

-0.0026 - 0.45231 

-0.0463 - 0.05841 

0.0095 - 0.05451 


-0.0163- 0.2401i 
0.0075 + 0.2518i 
0.2396 + 0.0827i 
0.2947 


0.0439 + 0.2918i 
-0.0356 - 0.2367i 
0.2487 + 0.0009* 
0.2211 


0.1131 + 0.0665* 
-0.0967 + 0.1000* 
-0.0398- 0.1196* 
0.1104 


-0.0341 - 0.0538* 
-0.0793 + 0.0196* 
0.0094 + 0.0408* 
0.0427 


-0.0233- 0.0012i 
-0.0148- 0.0912i 
0.0080 + 0.0213* 
0.0031 


-0.0223 + 0.0155* 
-0.0171 - 0.0985* 
0.0106 + 0.0049* 
0.0042 


-0.4516 - 0.0735* 
0.0336 + 0.0064i 
-0.0033 + 0.0768* 
0.4842 


0.4662 + 0.0508* 
0.0200 + 0.0379* 
0.0111 - 0.0464i 
0.5065 


-0.0026 + 0.4523* 
-0.0463 + 0.0584i 
0.0095 + 0.0545* 
0.4454 



PpA = 


Ppa = 


Pp.6 = 


Pp,7 — 


Pp.8 = 


Pp.9 = 


Pp.10 


Pp.ii = 


Pp.12 


Pp.13 


0.2807 

0.2569 + 0.0406* 

0.0456 - 0.2118* 

-0.0623 + 0.2607* 

0.2569 - 0.0406* 

0.2543 

0.0401 - 0.1770i 

-0.0304 + 0.2331* 

0.0456 + 0.2118* 

0.0401 + 0.1770i 

0.1256 

-0.2027 - 0.0279* 

-0.0623 - 0.2607* 

-0.0304 - 0.2331* 

-0.2027 + 0.0279* 

0.3394 

0.0260 

-0.0969 + 0.0186* 

0.0380 + 0.0730i 

0.0156 + 0.0018* ■ 

-0.0969- 0.0186* 

0.5035 

0.0389 - 0.4553* 

0.0051 - 0.0374i 

0.0380 - 0.0730i 

0.0389 + 0.4553* 

0.4676 

-0.0051 - 0.0057* 

0.0156- 0.0018* 

0.0051 + 0.0374i 

-0.0051 + 0.0057* 

0.0029 

0.0367 

0.0502 + 0.0923* 

0.0734 + 0.0820i 

-0.0044 + 0.0326*' 

0.0502 - 0.0923* 

0.4500 

0.4397 - 0.0152i 

-0.0288 - 0.0230* 

0.0734 - 0.0820* 

0.4397 + 0.0152i 

0.5064 

-0.0430 - 0.0023* 

-0.0044 - 0.0326* 

-0.0288 + 0.0230* 

-0.0430 + 0.0023* 

0.0069 

0.0062 

0.0557- 0.0422i 

0.0047 + 0.0036* 

-0.0086 - 0.0005*' 

0.0557+ 0.0422i 

0.9818 

0.0283 + 0.0662i 

-0.0441 -0.1031* 

0.0047 - 0.0036* 

0.0283 - 0.0662i 

0.0032 

-0.0049 + 0.0027* 

-0.0086 + 0.0005* 

-0.0441 + 0.1031* 

-0.0049 - 0.0027* 

0.0087 

0.5550 

-0.0622 + 0.0204i 

0.0308 - 0.0306* 

0.0026 - 0.4514i ' 

-0.0622 - 0.0204i 

0.0135 

-0.0451 + 0.0076* 

-0.0050 + 0.0533* 

0.0308 + 0.0306* 

-0.0451 - 0.0076* 

0.0043 

-0.0039 - 0.0220* 

0.0026 + 0.4514i 

-0.0050 - 0.0533* 

-0.0039 + 0.0220* 

0.4271 

0.0089 

-0.0219- 0.0731* 

-0.0027+ 0.0214i 

-0.0215 - 0.0387*' 

-0.0219 + 0.0731* 

0.5684 

-0.4398 - 0.0689* 

-0.0065 - 0.0551* 

-0.0027- 0.0214i 

-0.4398 + 0.0689i 

0.4209 

-0.0068 - 0.0044i 

-0.0215 + 0.0387* 

-0.0065 + 0.0551* 

-0.0068 + 0.0044i 

0.0018 

0.0058 

0.0406 - 0.0232i 

0.0029 + 0.0025* 

-0.0090 - 0.0064i' 

0.0406 + 0.0232i 

0.5133 

0.0500 + 0.4689* 

-0.0279 - 0.0088* 

0.0029 - 0.0025* 

0.0500 - 0.4689* 

0.4753 

0.0465 + 0.0408* 

-0.0090 + 0.0064i 

-0.0279 + 0.0088* 

0.0465 - 0.0408* 

0.0056 


0.2918 0.2395 - 0.0021i -0.0242 + 0.2258z 

0.2395 + 0.0021i 0.1992 -0.0291 + 0.2209z 

-0.0242 - 0.22581 -0.0291 - 0.2209z 0.2948 

0.0052 + 0.2556Z 0.0132 + 0.2382i -0.2469 + 0.0169z 


0.0052 - 0.2556i 
0.0132 - 0.2382i 
-0.2469 - 0.0169* 
0.2142 


0.1923 

-0.2541 + 0.0502i 
-0.0467 + 0.2470* 
0.0350 + 0.2186* 


-0.2541 - 0.0502i 
0.3002 

0.0458 - 0.2493* 
0.0075 - 0.2518* 


-0.0467 - 0.2470* 0.0350 - 0.2186* 
0.0458 + 0.2493* 0.0075 + 0.2518* 
0.2127 0.2396 + 0.0827* 

0.2396 - 0.0827* 0.2947 


0.3064 

-0.2359 + 0.0153* 
0.0152 - 0.2270* 
0.0170 - 0.2694i 


-0.2359 - 0.0153* 0.0152 + 0.2270i 
0.1884 0.0105 - 0.2216* 

0.0105 + 0.2216* 0.2840 

-0.0356 + 0.2367* 0.2487 - 0.0009* 


0.0170 + 0.2694i 
-0.0356 - 0.2367* 
0.2487 + 0.0009* 
0.2211 
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Ppm — 


0.3803 -0.0814 - 0.3887i 

-0.0814 + 0.38871 0.3748 

-0.1318-0.10511 -0.0803+ 0.147H 
0.1286- 0.0801Z -0.0967-O.lOOOz 


-0.1318+ 0.1051Z 0.1286+ 0.0801Z 
-0.0803- 0.1471Z -0.0967+ 0.1000i 
0.1345 -0.0398- 0.1196Z 

-0.0398+ 0.1196Z 0.1104 


0.5276 0.1680+ 0.4105Z 

0.1680- 0.4105Z 0.3819 

-0.0748 - 0.0527Z 0.0050 - 0.0826z 
-0.0633 + 0.0764i -0.0793 - 0.0196z 


-0.0748 + 0.0527Z -0.0633 - 0.0764i 
0.0050 + 0.0826Z -0.0793 + 0.0196z 
0.0478 0.0094 + 0.0408z 

0.0094 - 0.0408Z 0.0427 


0.0207 -0.0261 - 0.1296Z -0.0029 + 0.0401z -0.0181 - 0.0374i' 

-0.0261 + 0.1296Z 0.6312 -0.4367 - 0.0605z -0.0148 - 0.0912i 

-0.0029 - 0.0401Z -0.4367+ 0.0605Z 0.3450 0.0080 + 0.0213z 


-0.0181 + 0.0374i -0.0148+ 0.0912i 0.0080 - 0.0213z 0.0031 


Ppm = 


0.0202 -0.0269- 0.1385Z 

-0.0269+ 0.1385Z 0.7415 

-0.0033 - 0.0486Z -0.3803 + 0.0547z 
-0.0181 + 0.0334i -0.0171 + 0.0985Z 


-0.0033 + 0.0486Z -0.0181 - 0.0334i 
-0.3803 - 0.0547Z -0.0171 - 0.0985z 
0.2341 0.0106 + 0.0049Z 

0.0106 - 0.0049Z 0.0042 


Coefficient matrices 

All the analyzed tomographies are based on solving the linear-system problem 

Ax = b, 

where A is the coejficient matrix, b is the observation vector, and x = vec(p) is a real vector describing the unknown 
state p , i.e., 

x = vec{p) = [pii,Repi2,Impi2,Repi3,Im/9i3,...,p44]^. 

Thus, a two-qubit density matrix p is represented as a real vector x = (xi, ...,xie) with its elements given as follows 

XI X2 + iX3 X4 + ZXs Xq + ZXy 

X2 - iX3 Xs Xg + iXiQ Xu + ixi2 
X4. — ZX5 Xg — ZX3Q X33 X34 ZX35 

Xq — iX'i' Xii — iXi2 Xi4 — iX\Q '^IQ 

The coefficient matrices depend on the choice of the equations used for reconstructing a given density matrix. Below 
we list the transposed (for typographic reasons) coefficient matrices for the four analyzed tomographic protocols: 

■ 0000020020000000 ‘ 

00000 0 -2 002 0 0 0 0 0 0 

00020 0 0 000 -2 0 0 0 0 0 

0002000000200000 
00000 0 -2 00 -2 0 0 0 0 0 0 

00000 -2 0 020 0 0 0 0 0 0 

0000 -2 0 0 000 0 2 0 0 0 0 

. T _ 0000 -2 0 0 000 0 -2 0 0 0 0 

02000 0 0 000 0 0 0 -2 0 0 

00-20000000000020 
10000 0 0-100 0 0-10 0 1 

10000 0 0 100 0 0 -1 0 0 -1 

0200000000000200 
00 -2 00 0 0 000 0 0 0 0 -2 0 

10000 0 0 -1 00 0 0 1 0 0 -1 

1000000100001001 
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A 


T _ 

s — 


A^j = 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4 

0 

0 

0 

0 

0 

0 

0 

0 

4 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

0 

0 

0 
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Observation vectors 

The observation vectors correspond to photon coincidence counts. In reality we measure disturbed quantities 
b = b + S b instead of b. The observation vectors are column vectors. For convenience we arrange them in arrays, 
where each column corresponds to one of the 17 reconstructed states. 


Pi 

P2 

Pl7 

h.l 

§1,2 

blS7 

^2.1 

^2,2 

^2,17 

bN,i 

bN,2 

• &Af,17 


Note that the values of b listed below are not normalized and cannot be interpreted as probabilities. The elements of 
each vector b were registered over 5 seconds. This means that if an element of 6 is a sum or a difference of n projectors 
the measurement for each of the n projectors took 5/n seconds. In this way the measurements for observation vectors 
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of the same length take the same amount of time. To obtain the frequencies we can divide these values by the total 
number of photon coincidences counted or by a sum of coincidences counted for a set of projectors forming a basis. 
The set of such projectors is not unique. In our calculations we use the unnormalized coincidences and normalize the 
reconstructed density matrices. 


Standard 36 state tomography 
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The rows of the observation vector bs correspond to the following consecutive projectors: \HH){HH\, \HV){HV\, 
\HD){HD\, \HA){HA\, \HL){HL\, \HR){HR\,\VH){VH\,\VV){VVl \VD){VDl\VA){VAl\VL){VLl\VR){VRl 
\DH){DH\, \DV){DV\, \DD){DD\, \DA){DA\, \DL){DL\, \DR){DR\,\AH){AH\, \AV){AV\, \AD){AD\, \AA){AA\, 
\AL){AL\, \AR){AR\, \LH){LH\, \LV){LVl \LD){LD\, \LA){LAl \LL){LL\, \LR){LR\, \RH){RH\, \RV){RV\, 
\RD){RD\, \RA){RA\, \RL){RL\, \RR){RR\. 
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JKMW tomography 
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The rows of the observation vector bj correspond to the following consecutive projectors: \HH){HH\, \HV){HV\, 
\HD){HD\, \HL){HL\, \VH){VH\, \VV){VVl \VD){VDl \VL){VLl \RH){RHl \RV){RVl \RD){RDl \RL){RLl 
\DH){DH\, \DV){DV\, \DD){DD\, \DR){DR\. 


MUB-based tomography 
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1413 

1354 

1335 

1798 

1108 

1501 

950 

1197 

1398 

33 

1001 

978 

1447 

953 

1452 

1074 

31 

601 

1085 

718 

940 

1611 

907 

1238 

1253 

1422 

11 

1235 

1018 

1069 

1099 

756 

1067 

1769 

9 

852 

664 

1581 

1216 

1383 

1643 

1302 

1692 

1764 

1624 

776 

995 

1132 

778 

934 

21 

6227 

991 

1250 

1270 

1697 

1865 

1224 

1120 

722 

2102 

1407 

1433 

1210 

1150 

1200 

1057 

5667 

266 

2409 

2782 

97 

5284 

2720 

1181 

135 

134 

75 

2337 

16 

23 

1323 

958 

1438 

2600 

1867 

23 

21 

5418 

98 

2143 

874 

119 

27 

133 

2943 

133 

31 

1279 

1095 

1085 

957 

2362 

277 

268 

32 

55 

87 

1558 

3494 

4939 

1892 

88 

175 

2121 

1091 

868 

1002 

826 

1537 

586 

623 

157 

27 

109 

478 

1725 

253 

1887 

32 

4213 

2268 

1328 

775 

1148 

2240 

853 

4544 

4390 

2022 

1093 

184 

120 

182 

1636 

1011 

2320 

786 

1912 

4230 

57 

48 

607 

2921 

944 

945 

1245 

851 

111 

130 

41 

646 

637 

2232 

936 

2583 

52 

4028 

9 

693 

1795 

931 

998 

1404 

1117 

1788 

40 

2693 

1246 

990 

60 

1415 

71 

43 

70 

3875 

2676 

868 

1920 

1922 

1355 

1652 

3361 

4078 

1349 

1156 

613 

115 

862 

81 

13 

84 

54 

2240 

1613 

868 

921 


The rows of the observation vector 5 m correspond to the following consecutive projectors: \DH){DH\, \DV){DV\, 
\AH){AH\, \AV){AV\, \LD){LD\, \LA){LA\, \RD){RD\, \RA){RAl \VR){VR\, \VL){VL\, \HR){HR\, \HL){HL\, 
|$+)(<i>+|, |$-)($-|, |d/+)(vE-+|, \^-){^-\, ^{\DL)+z\AR)){{DL\-i{AR\), ^{\DL)-z\AR)){{DL\+z{AR\), U\DR) + 
i\AL)){{DR\ - i{AL\), ^{\DR) - i\AL)){{DR\ + i{AL\), where |4>±) = {\HH) ± \VV))/V2 and = {\HV) ± 
\VH))/V2. 
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Optimal tomography 


2727 

2575 

2844 

1448 

127 

193 

25 

2955 

40 

27 

1264 

809 

1231 

2762 

3831 

113 

112 

30 

58 

116 

1312 

2457 

2364 

3928 

72 

2555 

2375 

863 

1263 

757 

2722 

2773 

3452 

4102 

126 

43 

61 

648 

2282 

2660 

13 

23 

1892 

2199 

1277 

895 

1141 

977 

347 

1887 

1295 

2706 

2747 

2426 

1751 

14 

36 

35 

2274 

8 

26 

928 

1240 

888 

802 

310 

17 

23 

T08 

103 

444 

1325 

473 

263 

223 

334 

98 

188 

1037 

1069 

947 

594 

1219 

442 

449 

TM 

240 

241 

209 

94 

485 

169 

408 

328 

107 

9 

211 

61 

2823 

2980 

709 

766 

25 

104 

461 

235 

185 

385 

19 

164 

42 

13 

405 

496 

61 

957 

543 

46 

48 

231 

70 

131 

1092 

356 

434 

44 

163 

96 

44 

978 

1039 

942 

763 

383 

249 

269 

T8 

60 

52 

1045 

25 

226 

49 

20 

30 

215 

1069 

1008 

999 

289 

68 

43 

58 

429 

251 

297 

144 

28 

12 

44 

117 

20 

489 

73 

348 

3 

869 

296 

446 

27 

188 

108 

252 

457 

25 

454 

476 

26 

29 

429 

57 

31 

443 

702 

576 

84 

94 

36 

205 

348 

1202 

182 

121 

412 

284 

247 

40 

1032 

1059 

951 

726 

442 

498 

545 

62 

13 

TT 

540 

884 

2342 

2 

27 

2019 

73 

448 

46 

73 

706 

341 

1979 

1883 

1 

0 

14 

957 

2060 

393 

454 

7 

4 

2158 

978 

949 

799 

918 

634 

52 

56 

2660 

2593 

288 

153 

8 

53 

28 

302 

58 

4 

21 

68 

176 

821 

247 

427 

423 

48 

To 

2468 

nil 

38 

93 

44 

2501 

48 

42 

1153 

1010 

1172 

482 

390 

6 

85 


The negative elements of bo are marked with an overline. The rows of the observation vector bo correspond to the 
following consecutive measurements: \HH){HH\, \HV){HV\, \VV){VV\, \HD){HD\ — \HA){HA\, \HL){HL\ — 
\HR){HR\, \DH){DH\ - \AH){AH\, \LH){LH\ - \RH){RH\, \VD){VD\ - \VA){VAl \VL){VL\ - \VR){VR\, 
\DV){DV\-\AV){AV\, \LV){LV\-\RV){RVl |vl/+)(vl/+| _ |Mr-)(vl/-|, |T+)(T+| - |^-)(^-|, |$+)($+| - |$-)(cl>-1, 
|$+)($+| - |$-)($-|, where |$±) = {\HH) ± i\VV))/^/2 and |5'±) = {\HV) ± i\VH))/^/2. 


Pauli matrices based tomography 


1277 

1264 

132 

57 

133 

1144 

40 

443 

1037 

95 

52 

170 

55 

175 

242 

1244 

1102 

268 

497 

1241 

1129 

1115 

426 

134 

1222 

68 

1100 

1032 

984 

986 

825 

577 

63 

59 

82 

2 

357 

196 

80 

269 

98 

95 

9 

71 

84 

444 

102 

428 

17 

33 

38 

107 

107 

105 

39 

105 

117 

79 

69 

24 

58 

24 

83 

44 

830 

560 

49 

57 

143 

79 

1222 

246 

1107 

46 

432 

1182 

242 

1070 

75 

65 

96 

243 

23 

268 

244 

1247 

1264 

146 

264 

57 

1167 

74 

427 

940 

137 

74 

23 

43 

759 

248 

1145 

1002 

98 

138 

225 

1148 

269 

276 

244 

224 

472 

26 

1005 

1049 

932 

49 

420 

359 

407 

434 

68 

94 

55 

87 

455 

199 

64 

76 

15 

27 

40 

20 

745 

263 

140 

138 

45 

22 

196 

1186 

224 

245 

121 

455 

34 

44 

1054 

1039 

973 

454 

576 

93 

404 

137 

6 

269 

477 

60 

249 

90 

443 

154 

148 

32 

280 

33 

977 

1342 

413 

397 

1319 

1305 

1273 

310 

1150 

1199 

970 

1284 

1100 

1130 

13 

27 

55 

34 

255 

1302 

1316 

49 

39 

118 

90 

72 

35 

976 

183 

174 

44 

20 

46 

40 

926 

1487 

415 

724 

64 

82 

248 

140 

249 

19 

102 

476 

65 

202 

46 

34 

26 

440 

644 

50 

45 

292 

246 

28 

33 

32 

237 

79 

4 

174 

44 

41 

69 

29 

1846 

1638 

296 

370 

29 

47 

91 

242 

46 

113 

984 

158 

458 

44 

188 

200 

182 

54 

274 

367 

680 

1397 

1356 

1362 

1290 

1220 

1313 

1000 

1331 

1124 

1157 

1083 

1052 

1004 

1816 

1815 

1367 

1383 


The negative elements of bp are marked with an overline. The rows of the observation vector bp correspond to the 
following consecutive measurements: {\DD){DD\ + \AA){AA\) — {\DA){DA\ + \AD){AD\), {\DL){DL\ + \AR){AR\) — 
{\DR){DR\ + \AL){AL\), {\DH){DH\ + \AV){AV\) - {\DV){DV\ + \AH){AH\), {\DH){DH\ + \AH){AH\) - 
{\DV){DV\ + \AV){AV\), {\LD){LD\ + \RA){RA\) - {\LA){LA\ + \RD){RD\), i\LL){LL\ + \RR){RR\)- {\LR){LR\ + 
\RL){RL\), {\LH){LH\ + \RV){RV\)-{\LV){LV\ + \RH){RH\), {\LH){LH\ + \RH){RH\)-{\LV){LV\ + \RV){RV\), 
{\HD){HD\ + \VA){VA\) - {\HA){HA\ + \VD){VD\), {\HL){HL\ + \VR){VR\) - {\HR){HR\ + \VL){VL\), 

{\HH){HH\ + \VV){VV\) - {\HV){HV\ + \VH){VH\), {\HH){HH\ + \VH){VH\) - {\HV){HV\ + \VV){VV\), 

{\HD){HD\ - \HA){HA\) + {\VD){VD\ - \VA){VA\), {\HL){HL\ - \HR){HR\) + {\VL){VL\ - \VR){VR\), 

{\HH){HH\ - \HV){HV\) + {\VH){VH\ -\VV){VV\), {\HH){HH\ + \HV){HV\) + {\VH){VH\ + \VV){VV\). 





22 


Error analysis 


Estimated variances 


For all the tomographies the vectors of variances for the 17 measured states are given as matrices 


Pi P2 • • ■ Pn 

i>l r cr 2 ( 6 i,i) cr 2 (&i, 2 ) 0 - 2 ( 61 , 17 ) 

62 0 - 2 ( 62 , 1 ) 0 - 2 ( 62 , 2 ) 0 - 2 ( 62 , 17 ) 

: 

67 V _0-2(67v,i) 0 -2(6Ar,2) ... 0-2(6jv,l7) 


They can be approximated directly with 0-2(6) « 0-2(6) = 6 for all the tomographies except the optimal one and the 
tomography based on the Pauli matrices. For the optimal tomography the matrix of variances reads 


■ 2727 2575 2844 1448 127 193 25 2955 40 27 1264 809 1231 2762 3831 113 112 ' 

30 58 116 1312 2457 2364 3928 72 2555 2375 863 1263 757 2722 2773 3452 4102 

126 43 61 648 2282 2660 13 23 1892 2199 1277 895 1141 977 347 1887 1295 

2706 2747 2426 1751 14 36 35 2274 8 26 928 1240 888 802 310 17 23 

1352 1297 1445 1362 1299 1219 1936 1427 1250 1198 1088 1109 997 2842 3309 1900 2161 

1425 1456 1624 1433 1211 1207 1933 1515 1104 1102 1141 989 995 2844 3246 1700 2016 

1291 1199 1375 1015 1199 1395 23 1294 964 1158 1252 906 1190 1858 2177 955 669 

2 ,, . _ 1460 1427 1454 1184 1068 1203 16 1454 785 990 1058 1103 973 1936 2374 820 583 

O' ( 00 ) - ^421 1388 1196 1158 1076 1321 29 1154 981 1184 1104 1071 1035 935 284 930 660 

1369 1359 1204 1094 1225 1410 22 1118 998 1258 1026 1104 1070 900 305 968 691 

1419 1428 1233 1436 1124 1071 1770 1182 1236 1194 934 1236 834 1787 1563 1834 2080 

1295 1253 1196 1360 1352 1269 2108 1088 1328 1307 1055 1125 1037 1820 1674 1979 2329 

94 41 98 1018 2609 2596 1889 60 2194 2195 1209 821 1075 1533 1195 2565 2506 

44 53 146 1065 2112 2196 1772 86 2074 2276 1031 1114 862 1685 1628 2274 2316 

2757 2691 2431 1027 127 80 104 2640 74 27 1301 1026 1262 1778 2115 150 144 

2806 2583 2521 1249 160 164 105 2580 55 67 1183 1049 1227 1737 1964 306 274 


For the Pauli matrices based tomography the matrix of variances reads 


'1323 1313 1259 1206 1165 1201 903 1240 1093 1163 1069 1052 1054 1927 1799 1380 1377' 
1326 1329 1327 1185 1142 1257 853 1254 1022 1189 1075 1036 1032 1886 1877 1384 1340 

1355 1314 1304 1225 1162 1234 897 1238 1100 1176 1093 1072 1012 1823 1870 1395 1375 

1355 1314 1304 1225 1162 1234 897 1238 1100 1176 1093 1072 1012 1823 1870 1395 1375 

1392 1322 1323 1232 1228 1235 1089 1243 1089 1126 1042 1072 1016 1953 1952 1406 1455 

1370 1320 1322 1230 1225 1225 1032 1311 1051 1090 1009 1062 1000 1863 1965 1396 1407 

1378 1341 1325 1273 1210 1236 1062 1271 1057 1149 1056 1114 1005 1878 2024 1400 1456 

~ 1387 1343 1321 1261 1188 1270 983 1291 1116 1191 1097 1090 1016 1889 1797 1416 1411 

1397 1408 1414 1264 1218 1309 978 1317 1051 1180 1084 1047 1033 1872 1776 1334 1354 

1397 1356 1362 1290 1220 1313 1000 1331 1124 1157 1083 1052 1004 1816 1815 1367 1383 

1397 1356 1362 1290 1220 1313 1000 1331 1124 1157 1083 1052 1004 1816 1815 1367 1383 

1387 1343 1321 1261 1188 1270 983 1291 1116 1191 1097 1090 1016 1889 1797 1416 1411 

1397 1408 1414 1264 1218 1309 978 1317 1051 1180 1084 1047 1033 1872 1776 1334 1354 

1397 1356 1362 1290 1220 1313 1000 1331 1124 1157 1083 1052 1004 1816 1815 1367 1383 

1397 1356 1362 1290 1220 1313 1000 1331 1124 1157 1083 1052 1004 1816 1815 1367 1383 
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Estimated error radii 


For each tomography and reconstructed state we have estimated the maximum error R as described in the main 
text. Our results are summarized in the following matrix: 



Optimal 

MUB 

Pi 

■ 0.0983 

0.1475 

P2 

0.0997 

0.1554 

P3 

0.0987 

0.1506 

Pi 

0.0997 

0.1477 

Pb 

0.1051 

0.1622 

Pe 

0.0999 

0.1589 

P7 

0.1152 

0.2059 

P8 

0.0999 

0.1580 

II 

0.1083 

0.1712 

Pw 

0.1083 

0.1633 

Pll 

0.1130 

0.1526 

Pl2 

0.1135 

0.1597 

Pl3 

0.1179 

0.1570 

Pli 

0.0861 

0.1240 

Pw 

0.0864 

0.1281 

Pw 

0.0994 

0.1522 

Pl7 

0.0985 

0.1562 


Standard 

Pauli 

JKMW 

0.2183 

0.2407 

0.5712 

0.2221 

0.2448 

0.5357 

0.2231 

0.2422 

0.5297 

0.2130 

0.2304 

0.5044 

0.2344 

0.2566 

0.5767 

0.2291 

0.2449 

0.5553 

0.2637 

0.3222 

0.6065 

0.2299 

0.2455 

0.5850 

0.2507 

0.2682 

0.6135 

0.2388 

0.2668 

0.5807 

0.2281 

0.2511 

0.5829 

0.2293 

0.2563 

0.6009 

0.2318 

0.2631 

0.5850 

0.1509 

0.2044 

0.3907 

0.1587 

0.2128 

0.4412 

0.2171 

0.2479 

0.5441 

0.2168 

0.2531 

0.5229 


Note that the values are multiplied by a factor of 1.3 to compensate for underestimation of ||CT|j. Standard error is 
simply given by r = i?/2. 


Relative trace distances between the reconstructed states 

In order to compare the quality of the matrices reconstructed with different tomographic protocols, we have also 
calculated the relative trace distances for the respective states in each protocols. Here we omitted the JKMW protocol, 
because it provides the largest error radius. Having three relative distances (for the remaining three protocols) it is 
possible to visualize the relative distances between the matrices and their error radii on a plane. 


Pi 

P2 

P3 

Pi 

P5 

P& 

P7 

P8 

T = P9 
Pw 
Pll 
Pl2 
Pl3 
Pl4 
Pl5 
Pie 
Pl7 


T{po,Pm) 

0.1415 

0.1462 

0.1018 

0.1295 

0.0818 

0.1234 

0.0806 

0.1155 

0.1613 

0.1397 

0.0590 

0.1000 

0.0896 

0.0688 

0.0718 

0.1506 

0.1576 


T{po,Ps) 

0.1004 

0.0798 

0.1130 

0.1786 

0.1684 

0.1176 

0.0483 

0.1150 

0.1245 

0.0956 

0.0515 

0.1074 

0.1000 

0.0425 

0.0790 

0.1270 

0.1278 


T{pm,Ps) 

0.1203 

0.1048 

0.1362 

0.1967 

0.1924 

0.0818 

0.0990 

0.1541 

0.0998 

0.1404 

0.0591 

0.0960 

0.0998 

0.0686 

0.0813 

0.1311 

0.1179 
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FIG. 4. Relative distances between points representing the reconstructed density matrices and their corresponding circles of 
the maximum errors R optimal (O), standard (S), MUB-based (M) tomographies. The 12 states pn are given in the units 
of trace distance. The states can be approximated by using Eq. 0 with Pn = \ilJn}{4’n\. The absolute positions of the three 
points are irrelevant. All the states for n = 1,12 are fully entangled except pj. The states of for n = 13,17 are partially 
entangled or separable (n = 13, 14). The ideally reconstructed state lies in the intersection of the error circles of radius R. 





















